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Cosmological solutions of Einstein's equations for equilibrium statistical systems of particles 
with scalar interaction are investigated. It is shown that the scalar field can effectively change the 
state equation of a statistical system, that leads to the possibility of secondary acceleration of the 
cosmological expansion. 
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I. INTRODUCTION 

To the moment it has been published many papers devoted to the secondary acceleration of the Universe. In 
order to solve the problem of the secondary acceleration of the Universe it is often proposed in many papers to 
sufficiently change the fundamental principles of physics. However, recently there are indications that complicated 
multi-component, classical physically systems can also bring to the secondary acceleration of the Universej^. In this 
case there is no need to revise the fundamental principles of physics. Some indications of the possibility of such 
behavior of the systems with scalar ineraction were given in the paper Scalar fields were introduced in general 
relativistic statistics by one of the Authors in the early 1980 years 0|, Q. In these papers, based on the kinetic theory, 
a system of particles with scalar interaction is obtained. 



II. MACROSCOPIC DENSITIES 



The full degeneration condition supposes: 



^ ^ oo. 



(1) 



(/i - chemical potentials, 9 - temperature) local equilibrium distribution function of Fermi-system is a step function 
(see for example Q). In this case, the local equilibrium distribution function has the form 



where X+(z) is a Heaviside function; 



m* = m 



9*1, 



(2) 



(3) 



are the particles effective masses, g is a scalar charge of the fermions, $ is a scalar field potential, m is a vacuum mass 
of the fermions. Therefore the integration of macroscopic densities are representable in elementary functions [3]: 



£/ = ^[V'x/r+v^(i+2^2) 



pf 
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[7/;Vl + V'2(2V''-3) 



+31n(?A + v/T+V^)]; 
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ef+Pf = ^i^'Vl + i^; (7) 



a = T/ 



q ■ ml 
27r2 

where 



^^Y+^ -M'<p + ^/i + i?)]- (8) 



V'^p.p/lm*! (9) 

is a ratio of the Fermi momentum to effective mass, ef,pf,a is an energy density, pressure and scalar charge density 
of fermions, respectively In this case the self-consistent equation of massive scalar field with a mass of fermions /i 
i: 

+ = -Aira. (10) 
The fermion number density is connected with the Fermi momentum by relation 

nix) = ^ =^ Pf = {SttMx))'^ ■ (11) 

Thus, the variable tp can be expressed in terms of two scalars - the particles number density, nix), in the own frame 
of reference and the scalar potential, 

i^;rM^_ (12) 

III. SELF-CONSISTENT SYSTEM OF EQUATIONS FOR A FLAT FRIEDMANN MODEL 

Let's consider a cosmological situation where matter is represented by only a degenerate Fermi system with scalar 
interaction of particles and the associated scalar field, described by the equation ([T0| (see fl'|). As the metric we choose 
the Spatially Flat Friedmann Model [5|]): 

ds^ =dt^ -a^it)idx^ +dy^ +dz^), (13) 
in which the independent Einstein equations have the form: 

3^ = Stte; (14) 

3^ = -^. (15) 
a e + P ^ ' 

In this metric the conservation law for particles takes form: Q: 

dt\/~Hn = 0, 

from where, taking into account epression pip we obtain a momentum integral instead of the energy integral fll]: 

apF = Const. (16) 

Supposing further $ = $(t), we obtain e = eit) p — pit) and the structure of the total scalar field EMT in the form 
of the EMT of a perfect fluid with the macroscopic velocity = 5\, energy density es(i) and pressure Psit): 

ss^^i^' + f^l'i''); Ps = ^il^'-f^l'^'). (17) 
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Differentiating (fT7|) and substituting tlie resuits in it can be seen, tliat tliis equation is a consequence of field 
equations ([TU)) . which in metric takes form: 

Thus, for an independent equations one can chose the field equation ^5)1 and one of Einstein's equations p^ . which 
with the condition (d > 0) can be written as: 



a 8n , , 

a - V T^' ^''^ 

where 

e = £/ + 

is a total energy density of fermions and a scalar field. 

Let us carry out an analysis of possible cosmological evolution of the degenerate scalar charged fermions system. 
From the expressions for the macroscopic densities ([4]), ([5]) we obtain: 

Pf = l^f 
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4V'Vr+V^-31n(V'+ Vl + V^) • (20) 



We can show, that the expression in square brackets on the right side (|20l) is a nonnegative for nonnegative "0^0- 
Therefore the strict inequality is satisfied: 

0<Pf<lef. (21) 
Similarly from p7|) . we find the relation for the scalar field: 

Ps = - ^A*'*', (22) 
Also for the scalar field the strict inequality is satisfied: 

Ps < ^e.. (23) 

Now we find coefficient of barotropic, as a coefficient k in the linear relation between total pressure and energy density 
of matter: 

p = K£. (24) 

As a result we find: 
where: 

G(V') = A^P^/lT^ ~ 3 InitP + VTT^). 
So, we can show, that for the concerned system the following relation is always fulfilled: 

-1 < ^{t) < 1 (26) 

Thus, for the cosmological acceleration Universe, ft, that consists of degenerate fermions and scalar field, 

ad 1 
2 

we obtain the restriction: 



n='^ = --{! + 3k), (27) 



i<n<i. (28) 
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IV. NUMERICAL SOLUTION OF EINSTEIN'S EQUATIONS FOR A DEGENERATE PLASMA WITH 

SCALAR FIELD 



Attempts of direct numerical integration of Einstein and Klein-Gordon equations (|14p . (|10p in most cases don't give 
results owing to nonlinear character of these equations and ambiguity of radical functions and logarithms in the right 
side of equations. Therefore, for the numerical solution of Einstein's equations the expressions for the macroscopic 
densities were extrapolated by means of elementary functions. As it can be seen from (IH - ([8|), all these expressions, 
taking into account the integral (jl6p up to multiplication by a conformal factor 1/a*, - 



N 4 

aoPp 



atp 



(29) 



are elementary functions of a single dimensionless function ip. This fact allows us to find the interpolation expression 
for the corresponding conformal densities: 



Pf = a Pf, if 



4 

a a. 



(30) 



Interpolation of analytic functions for the energy density, pressure and charge density on the range under consideration 
can be written as: 

Pf = ^^e-^H + (31) 
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-24)-2tlj- 



)( 



15 

38 I _±_ 

15 3-02 



(32) 
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15 



~2i/' 



2^-2 



(33) 



1 + ■ 

It is obvious that due to the large number of essential parameters of the model, (4 parameters, Pf,m, ^,q)^ and 
initial conditions, (2 independent conditions, <i>(0), $(0)), a considered cosmological model is extreme rich in types of 
behavior. Therefore let us consider the main ones. Here are the results of numerical solutions of Einstein and Klein - 
Gordon equations for the degenerated Fermi system in math package Mathematica v7. 
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Fig. 1. Evolution of the barotropic coefficient, depending on 
mass of scalar bosons. Frequent dotted line: /i = 0,3; dotted line: 
/Lt = 0, 1; rare dotted line - = 0, 2; normal line: /i = 0, 25; thick 
line: = 0, 35. Everywhere: pp = 0, 01; to = 1; q = 0, 3; $(0) = 
1; $(0) = 0. 
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Fig. 2. Evolution of the cosmological acceleration, il, depending 
on mass of scalar bosons. Frequent dotted line: fj, = 0,3; dotted 
line: = 0, 1; rare dotted line: /i = 0, 2; normal line: fi = 0,25; 
thick line: /i — 0,35. Everywhere: pp = 0,01; m = 1; q = 
0,3;«'(0) = 1; $(0) = 0. 




Fig. 3. Evolution of the cosmological acceleration, fl, depending 
on mass of scalar bosons. Frequent dotted line: /i = 1; dotted 
line: /i = 1,2; normal line: fi — 2; thick line: fj, = S. Everywhere: 
PF = 0,01; m = l; q = 0, 3; $(0) 1; $(0) = 0. 

Let us present the solutions which transgress to the accelerated cosmological phase with maximum acceleration. These 
solutions correspond to very small values of the fundamental charge. 



Fig. 4. Evolution of the cosmological acceleration, il, depending 
on mass of fermions. Frequent dotted line: to = 0; dotted line: 
TO = 0, 001; rare dotted line: to = 0, 01; very rare dotted line: to ~ 
0.03; thin line: to = 0, 05; normal line: to — 0, 1, thick line: to — 1. 
Everywhere: = 0,01; /i = 0,3; g=l;$(0) = l; $(0) = 0. 
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Fig. 5. Evolution of the cosmological acceleration, where "pp = 
0; TO=1; /i = 0,001; g ^ 0. 
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Fig. 6. Evolution of tlie cosmological acceleration, fl where pp — 
0; m = 1; /i = 0, 1; q ^ 0. 

V. CONCLUSION 

Thus, we can resume: 

1. Evolution of the cosmological model depends heavily on the fundamental constants and m has many types of 
behavior; 

2. There are ranges of fundamental constants and initial conditions, in which the cosmological models come out later 
on stable acceleration; 

3. Access to the phase later acceleration can be smooth or accompanying by the acceleration oscillations; 

4. Final acceleration is constant and can vary in the range —1 < 57 < 1. 
Therefore, constant acceleration (deceleration) is of a general nature. 
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